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Translation is an important process for prokaryotic and eukaryotic cells to produce necessary 
proteins for cell growth. Numerious experiments have been performed to explore the translational 
properties. Diverse models have also been developed to determine the biochemical mechanism of 
translation. However, to simplify the majority of the existing models, the frameshifting of ribosome 
along the mRNA transcript is neglected, which actually occurs in real cells and has been extensively 
experimentally studied. The frameshifting of ribosome evidently influences the efficiency and speed 
of translation, considering that the peptide chains synthesized by shifted ribosomes will not fold 
into functional proteins and will degrade rapidly. In this study, a theoretical model is presented 
to describe the translational process based on the model for totally asymmetric simple exclusion 
process. In this model, the frameshifting of the ribosome along the mRNA transcript and the 
attachment/detachment of the ribosome to/from the main body of mRNA codons during translation 
elongation process, are explicitly included. The results show that, with ribosome frameshifing, the 
speed of correctly synthesized peptide chains may increase first and then decrease with both the 
translation initiation rate a and the ribosome detachment rate aid. This results indicates that 
regulating the translation process to reach maximal synthesized speed of proteins is theoretically 
feasible. Traffic-related problems of ribosome motion along the mRNA transcript are also addressed 
theoretically. Depending on parameter values, shock wave (or domain wall) may exist for ribosome 
probabilities along the mRNA. 


I. INTRODUCTION 

In living cells, proteins are synthesized by ribosome through the translation process with the information coded in 
the template messenger RNA (mRNA). The mRNA is composed of a sequence of codons, and each codon consists of 
three nucleotides. Given the four kinds of nucleotides in mRNAs, a total of 4^ = 64 kinds of codons exits. Except 
for the stop codon, which is used to finish the translation process, each kind of codon corresponds to a certain kind 
of amino acid, although each kind of amino acid may correspond to several kinds of codons [1]. During translation, 
an amino acid is transferred by the transfer RNA (tRNA). With the binding of tRNA to mRNA through one of its 
end, which has complementary anticodon sequence to that of the mRNA, the amino acid carried by the other end 
of tRNA is then chained together into a polypeptide as the mRNA passes through. A schematic description of this 
process is found in Fig. 1. 
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Translation is followed by molecular machine ribosome, which includes a three step process, that is, initiation, 
elongation, and termination [2]. The 308 subunit of the ribosome contains three tRNA binding sites, which are 
denoted by A, P, and E (Fig. 1). In the initial step, the first tRNA and the start codon of mRNA will form a 
codon-anticodon duplex in the P site of the 308 subunit, with the 8hine-Dalgarno (8D) sequence on the 5’ end of 
the mRNA binding with the 3’ end of the 168 rRNA (an anti-8D sequence in the 308 subunit of a ribosome) [3-7]. 
Meanwhile, the 508 subunit of ribosome will attach to the mRNA from the other side. During translation, mRNA 
will be kept in the channel between these two subunits. During elongation, each binding site of the 308 subunit of 
ribosome will be bound by one mRNA codon, while the tRNA anticodon end will bind to the mRNA codon in the A 
site. When the ribosome steps forward along the mRNA, the codon from the A site, together with the tRNA attached 
to it, will move to the P site. Then, the amino acid carried by the tRNA will be added to the tail of the synthesizing 
peptide chain in the P site. Finally, when the complex of mRNA codon and tRNA moves to the E site, the tRNA 
will detach from the mRNA codon and the peptide chain (Fig. 1). No tRNA that contains an anticodon that is 
complementary to the stop codon exists. In the termination process, the stop codon will be bound by one release 
factor, which helps the release of the peptide chain and ribosome from the mRNA, which completes the translation 
process. 

The translation process described above is only applicable for ideal cases, in which the ribosome always moves 
forward codon by codon, and it always reads three nucleic acids from the same codon in each elongation step. 
However, in actual cells, the ribosome may shift to one nucleic acid upstream or downstream along the template 
mRNA in a low frequency [8]. This phenomenon is called frameshifting. Two mRNA codons are considered to be in 
the same frame if and only if the number of nucleic acids between theses codons is divisible by 3. In the following, 
the mRNA start codon is called a 0 frame codon (or a codon in 0 frame). Thus, all codons that are in the same frame 
with the start codon are also 0 frame codons (or codons in 0 frame). A codon that is one nucleic acid upstream to 
a 0 frame codon is called a —1 frame codon (or a codon in —1 frame). Meanwhile, a codon that is one nucleic acid 
downstream to a 0 frame codon is called a -1-1 frame codon (or a codon in -fl frame). Considering that every codon 
consists of three nucleic acids, any mRNA codon must belong to one of the three frames, i.e 0, —1, and -1-1 frames. 
If one mRNA codon read by ribosome belongs to the 0 frame (or ±1 frame), then the mRNA is said to be read in 
the 0 frame (or ±1 frame). Therefore, each mRNA has three reading frames. With frameshifting of ribosome, the 
reading frames of mRNAs are distinct. To obatin intuitive impressions of frameshifting, one schematic diagram of —1 
frameshifting is depicted in Fig. 1. 

Translation is generally initialized from the start codon, so, in the initiation of translation, mRNA is usually in the 
0 frame. During the elongation period, frameshifting may occur at some codons, so the mRNA will change into the 
— 1 frame or -|-1 frame. Regardless of the reading frame, translation will continue until the ribosome meets the stop 
codon. However, peptide chains produced by ribosome in ±1 reading frames are different from those produced by 
correct translations in the 0 reading frame, and these chains are usually nonfunctional, and will degrade rapidly. 

Experimental findings have shown that the frequency of frameshifting is only 10“^ to 10“^ per codon. However, 
the frequency of —1 frameshifting increases for a slippery sequence with the motif X XXY YYZ [9]. The 3’ secondary 
structure and the 5’ internal 8D sequence in the mRNA can further increase the —1 frameshift frequency [10, 11]. 
Numerous studies have been performed to describe the property and mechanism of frameshifting [12-20]. A general 
mechanistic and conformational framework for —1 frameshifting is recently presented in [9] to attempt to elucidate 
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the mechanism of frameshifting. Structural insights into +1 frameshifting are discussed in [21], and some possible 
mechanisms of frameshifting can also be found in [22]. 

Biophysically, translation can be roughly regarded as a totally asymmetric simple exclusion process (TASEP) see 
[23-26]. TASEP is one statistical physics model to describe the unidirectional hopping process of general particles 
along a one-dimensional lattice. In this model, the particle at site i of the main body of lattice hops forward to 
site i -I- 1 provided that site i -I- 1 is not occupied. The particle at the last site of lattice leaves the lattice into the 
environment at a given rate constant, and particles in the environment enter the first site of lattice provided that 
this site is unoccupied. If one mRNA with A^ -|- I codons is considered as a one-dimensional lattice of length TV -|- I, 
then the initiation process of translation corresponds to the entrance of the ribosome into the first lattice site. The 
elongation process of translation corresponds to the forward hopping process of ribosome along the main body of 
lattice, and finally the termination of translation corresponds to the leaving of ribosome from the last lattice site. 
Therefore, to a certain extent, TASEP is a reasonable model to describe the translation process in gene expression. 

Different types of TASEP models have been conceptualized in various kinds to describe corresponding biophysical 
or biochemical processes. In these models, particles may be allowed to detach from or attach to any site of lattice, 
particles may have multiple internal states, there may be different types of particles which hop along the same 
track, particles may be allowed to hop to the nonadjacent lattice sites, and particles may hop along multiple parallel 
lattices [27-36]. Nevertheless, no theoretical model to date has been presented to study the translational process with 
occasional frameshifting of ribosome along mRNA. Frameshifting influences on both the speed of correct translation 
(i.e., the product speed of correct peptide chain) and the accuracy of translation. Therefore, to determine the detailed 
translational properties, and explore the mechanism to regulate gene expression, a more reasonable model that includes 
frameshifting should be presented. In this study the translation process, which includes initiation, elongation and 
termination, will also be regarded as one TASEP of particle hopping along a one-dimensional lattice. Contrary to the 
usual TASEP, particles (i.e., ribosomes) in the translation elongation period may make one frameshifting stochastically. 
The corresponding peptide chains produced by shifted ribosomes will be nonfunctional and degrade quickly. Given 
the difference of ribosomes in correct and incorrect translation states, the proposed model is similar to the TASEP 
with two particle species. However, the proposed model is different because, during forward hopping, particles may 
change from one species to the other. Additionally, due to frameshifting, ribosomes in incorrect translation state do 
not stop at the supposed end of lattice. This is because that for ribosomes that translocate in the wrong mRNA 
frames (i.e., ±1 frame), the stop codon cannot be recognized as usual. This study was primarily designed to discuss 
the translational properties of by the modified TASEP model, which includes the frameshifting of ribosome, especially 
the dependence on model parameters of correct translation speed and the ratio of correctly produced peptide chains 
in all peptide chains. 

This study is organized as follows. The theoretical model will be presented in the next section, and then the 
numerical results on the translational properties with ribosome frameshifting will be provided in Section III. Finally, 
concluding remarks will be presented in the last section. 
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II. MODIFIED TASEP WITH RIBOSOME FRAMESHIFTING 

In the usual TASEP without ribosome frameshifting and attachment/detachment to/from the main body of mRNA, 
the model parameters that influence translation speed are as follows: (i) entrance rate of particles from the environ¬ 
ment to the first lattice site, i.e., the rate of ribosome binding to the mRNA start codon. This value depends on 
the concentration of free ribosomes, and the sequence of ribosome binding site (RBS); (ii) length of lattice (mRNA), 
which is usually denoted by N; (iii) leaving rate (3 of particles (ribosomes) from the last lattice site N (the stop 
codon) to environment; and (iv) forward hopping rate kE of particles in the main body of lattice, i.e., the elongation 
rate of translation. For simplicity, the particles are usually assumed to be able to hop forward along the track, and 
their forward hopping rates at any lattice site i are assumed to be the same. In this study, the usual TASEP had 
been modified to include frameshifting of ribosome as shown in the schematic depiction in Fig. 2. 

Given that the frequency of frameshifting is usually low, this study assumes that throughout the translocation 
along mRNA, each ribosome can only frameshift at most once. This assumption implies that only ribosomes in the 0 
frame can shift one nucleotide forward to 1 frame or backward to —1 frame. General models without this restriction 
can be presented similarly, but those models will complicate the following analysis. Actually, with more possible 
frameshitings or additional ribosome frameshiftings in the ±1 frame, the results obtained in this study will not change 
essentially. In this study, the codon with one nucleotide upstream or downstream of codon i is denoted by i— or i+. 
The rates of frameshifting from codon i to codon i— and codon i+ are denoted by fc_ and k+, respectively. 

Gonsidering that ribosomes only have one chance of framshifting, the ribosomes in ±1 frame will not shift back 
into the 0 frame. Thus, no stop codon will generally exist in the ±1 frame, and ribosomes in the ±1 frame will come 
across the stop codon and continue moving forward to the 3’ end of mRNA. This study denotes the number of codons 
between the usual stop codon and the 3’ end of mRNA by M, and they are numbered by (TV -|- 1)±, • • • , (TV -|- M)±. 
Thus, the length of track for ribosomes translocating in ±1 frame is TV -|- M, which is M codons longer than that for 
ribosomes translocating in the 0 frame (Fig. 2). The hopping rate of particles (elongation rate of translation) from 
codon T± to codon {i -|- 1)± are denoted by k^ and the probabilities of finding one ribosome at codon T± are denoted 

by qf- 

In addition to start and stop codons, this study also allows the attachment/detachment of ribosomes to/from any 
mRNA codon with rates denoted by Wa/w^. Attachment and detachment of ribosomes are also allowed at codon T± 
for any 0 < i < N + M with rates denoted by uj^ and respectively. Therefore, ribosomes at codon i for 1 < T < TV 
can be divided into two classes, namely, ribosomes with correctly and incorrectly synthesized peptide chains. If a 
ribosome reaches codon i from the start codon 0 without frameshifting, then the synthesized peptide chain will be 
correct. On the contrary, if a ribosome binds to mRNA at codon j for 1 < j < i, then the synthesized peptide chain 
will be incorrect. For 0 < i < TV, the probability of finding a ribosome at codon i with one correct or incorrect peptide 
chain is denoted by pi or qi. This study assumes that only ribosomes in the 0 frame with a correct peptide chain 
can frameshift, and the probability qo is always equals to zero. Considering that frameshifting of ribosomes with 
incorrect peptide chain will not change the correct translation speed essentially. Without loss of generality, except 
for the normal initiation and termination, this study assumes the absence of ribosome detachment/attachment at the 
start/stop codon. 

For simplicity, this study assumes that kE = k^,uja = and fc_|_ = =: fcg. This condition indi- 
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cates that for a ribosome regardless of frames, either in the 0 or in ±1 frame, its forward translocation rate and 
detachment/attachment rates are the same. Meanwhile, the rates of upstream and downstream frameshifting are also 
equal. 

At any time, only one of the three codons i, and i— can be occupied by ribosome. Thus, if codon i is occupied 
by one ribosome, then the codons *± will be empty. This phenomenon indicates that the ribosome at codon i always 
has the possibility to frameshift and then translocate to either codon i— or codon i+. For convenience, Wi denotes 
the probability that position i is unoccupied. One can show that (see Fig. 2), 

w = [ ^ “ 9+ - g* - 9” -Pz, 0 < z < Af, 

\l-qt-qi-q-, N + l<i<N + M. 

The probabilities pi of finding ribosome with correctly synthesized peptide chain at codon i are governed by the 
following equations: 


dpo/dt = -2ksPo + ctWo - ksPoWi, 

dpi/dt =-2ksPi +kEPi-iWi - ksPiWi+i - uidPi, l<i<N-l, (2) 

dpN/dt = -2ksPN + ksPN-iWN - Ppn- 


Similarly, the equations for probabilities q^ of finding ribosome at codon z± are as follows: 


dq^/dt = ksPo - kEqoWi + uJaWo - 

dqf / dt = ksPi + kEqt_iWi - kEqfW^+i + ojaWi - ojdqt, (3) 

dqf /dt = kEqf_iWi - kEqfW^+i + UJaW^ - OJdqf, N + l <i < N + M. 

Where Wat+m+i = 1- If M 0 then the probabilities qi of finding ribosome at codon i but with incorrectly synthesized 
peptide chain satisfy 


dqi/dt = kEqi-iWi - kEqiWi+i + UaWi - uJdq^ iov I <i <N + M and i ^ N,N + I, 
dqN/dt = kEqN-iWN — PqN, (4) 

dqN+i/dt = —kEqN+iWN+2 +^aWN+i — LOdqN+i, 

and similar equations can be easily obtained for M = 0. With the values of probabilities Pi,qi and q^, the effective 
rate of translation initiation can be obtained by cxeff '■= ct{^ — Po — Qo — Qq — %), and the effective translation 
rate, i.e. the rate of synthesizing correct peptide chains, can be obtained by Peff '■= Ppn- The ratio of Peff to 
Qfe//, r := Peff/oieff, IS One index to describe the effectiveness of translation along given mRNA template. Without 
ribosome frameshifting, attachment and detachment, the steady state value of r is equal to 1. For such simple cases, 
each ribosome that successfully begins its translation process from the start codon will finally complete its translation 
at the stop codon with one correctly synthesized peptide chain. 


III. RESULTS 

To illustrate the modified TASEP used in this study for the description of ribosome translocation along mRNA with 
frameshifting, typical examples of the probabilities Pi,qi, qf and their summation Pi = Pi + qi + qi^ + qf are plotted in 
Fig. 3. Where Pi = 1 — Wi is the probability that site i is occupied by one ribosome. These examples indicate that. 
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from the viewpoint of total probability Pi of ribosome, the modified TASEP can be regarded as a combination of 
two usual TASEPs for ribosome translocation along mRNA but without frameshifting. Where one involves ribosome 
translocation from lattice site 0 to lattice site N, and the other involves for ribosome translocation from lattice site 
N to lattice site N + M. 

Eor TASEP without ribosome frameshifting, attachment or detachment, i.e., kg = uia = oJd = 0, the theoretical 
studies in [37, 38] show a total of three possible phases, namely, low ribosome density phase, high ribosome density 
phase, and maximal current phase. The plots in Figs. 3(a-c) show that, such three phases also exist between lattice 
site 0 and lattice site N in the modified TASEP. An example for maximal current cases is found in Fig. 3(a), in which 
the total probability Pi is almost equals to 1/2 except for the sharp changes near sites 0 and N. The maximal current 
phase occurs when both a and /3 are bigger than 1/2. In TASEP, the current is defined by J = ksPiil — Pi)- The plots 
in Figs. 3(b) and 3(c) correspond to low and high density phases, respectively, with only sharp change at site N or site 
0. For these special cases, considering the absence of ribosome frameshifting, attachment or detachment, probabilities 
qi,qf are equal to zero, and pi = Pi- The current j = KEPi{i — Pi) of ribosome with correctly synthesized peptide 
chain is conversed along the mRNA region between the start and stop codons. Thus the effectiveness r of translation 
is equal to 1, and actually aeff and feff are both equal to the current j = KEPi{l — Pi) = kEPi{l — Pi) = J- 

For general TASEP, i.e., uja-,0Jd 0, previous studies have shown that domain wall (or shock wave) and sharp 
changes at both or only one of the boundaries may exist [24, 25]. The examples plotted in Figs. 3(d-i) show that, 
with frameshifting of ribosome, in the two regions [0, A^] and [N, N + M], domain wall and sharp change of probability 
at one or both of the two boundaries may be found. Meanwhile, the examples plotted in Fig. 3 show that, similar 
to the total probability Pi, three different phases may also exist for probability pi (for the cases uJa = uJd = 0, see 
Figs. 3(a-c)), domain wall and sharp change at boundaries may also exist (for general cases with nonzero attachment 
and detachment rates, see Figs. 3(d-i)). In all these figures, the total probability Pi decreases suddenly at codon 
Af + 1 = 151. This phenomenon is due to the existence of the stop codon at site N = 150. Instead of moving 
forward to codon A^ + 1 = 151, ribosomes at the stop codon will leave the mRNA template, and consequently the 
total probability Pi at codon TV + 1 = 151 will decrease rapidly. 

For translation with ribosome frameshifting, the main properties that should be determined are the influences of 
frameshifting on the rate of translation initiation aeff and the rate of synthesizing correct peptide chain Peff, as well 
as their ratio r. Generally, shifted ribosomes, i.e., ribosomes in ±1 frame of mRNA, will not leave the mRNA from the 
stop codon. The probability of establishing a ribosome at one codon, even if the codon is upstream of the stop codon 
(i.e., its codon number is less than TV), may be influenced by the ribosome probabilities between the stop codon and 
the 3’ end of the mRNA. Meanwhile, previous studies about usual TASEP show that ribosome probabilities between 
sites TV and TV + TIT depend on the length M of the untranslated mRNA region [26]. Thus, the ribosome probability 
at any mRNA codon may also depend on parameter TIT, and consequently, any properties of the modified TASEP 
model for translation with ribosome frameshifting obtained from the probabilities Pi,qi, and qf^ may be parameter 
Af-dependent. Considering the lack of experimental measured values for parameter TIT, one strategy used in this study 
is the selection of a sufficiently large value of TIT, such that all the interested properties of translation, especially the 
effective initiation rate aeff and termination rate Peff, reach their steady values. Or in other words, aeff and /3e// 
will not change with further increase in parameter TIT, i.e., the length TIT of the 3’ untranslated mRNA region is large 
enough to have as little influence as possible to the ribosome translocation in the translated region. The numerical 
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calculations show that, for the parameter values used in this study, the selection of one value of M that is larger than 
20 is sufficient (Fig. SI in [39]). In Fig. 3, M = 50 is used, and in Figs. 4-6, M = 25 is used. 

Without frameshifting, attachment, or detachment, ribosomes initiate their translations from the mRNA start 
codon and terminate their translations at the stop codon with correctly synthesized peptide chains. On the contrary, 
with frameshifting, attachment, and detachment, ribosomes may begin/end their translations at any mRNA codon. 
However, only ribosomes that initiate their translations from the start codon and leave the mRNA from the stop codon 
synthesize correct peptide chains. Thus, the synthesizing rate of correct peptide chain, which is called the effective 
translation rate /3e// := Ppn: will be less than the effective translation initiation rate Og// ■= o.{^—Po — Qo — Qq —Qq)- 
Their ratio r := j3eff/(y.eff will be one important biophysical index to describe the effectiveness of gene translation. 

In Fig. 4, the Cjlfectiwe translation initiation rate Og//, elective translation speed/?g//, and their ratio r = j3eff/cieff 
are plotted with the change in translation initiation rate a. The plots in Figs. 4(a,b) show that, without frameshifting 
ks = 0, both ttg// and /3g// increase with a. By contrast, for nonzero attachment rate uia, rates Og// and /3g// will 
reach their up limits rapidly. Given that the new ribosome attachment will quicken the saturation of the ribosome on 
the mRNA. Fig. 4(c) shows that, if ivd = 0, then the ratio r = fif,ff/aeff is equal to I. In such cases, all ribosomes 
that begin their translations from the start codon will finally end their translations at the stop codon. Generally, 
r < 1 holds, and r decreases with initiation rate a and finally tends to one low limit value. For kg > 0, i.e., translation 
with ribosome frameshifting, the plots in Fig. 4(d) indicate that the effective translation initiation rate Og// increases 
also with Of, while the plots in Fig. 4(e) show that the effective translation rate /3g// increases first and then decreases 
with rate a. The reason that /3g// decreases with large values of a is that for large a, the probability of finding a 
ribosome with incorrect peptide chain will increase, so the translocation speed of ribosome with correct peptide chain 
will decrease because of the block of ribosome with incorrect peptide chain. Finally, the plots in Fig. 4(f) indicate 
that for ks > 0, the ratio r will always be less than I and decreases with translation initiation rate a. The plots in Fig. 
4 show that the main differences that resulted from ribosome frameshifting are as follows: the effective translation rate 
Peff may decrease with translation initiation rate a (see Fig. 4(b,e)), and the effectiveness r of translation decreases 
much rapidly with rate a (see Fig. 4(c,f)). 

Meanwhile, the plots in Figs. 5(a-c) show that cteff, Peff-, and their ratio r increase with the translation termination 
rate /3, and tend to corresponding limit constants with large /3. These results imply that, with high leaving rate fi 
of ribosome from the stop codon, more correct peptide chains is synthesized, and the effective initiation rate also 
increases. Given the increase in leaving rate (3 but hxed rates of frameshifting and detachment, ribosomes with 
correct peptide chain rapidly translocate along the mRNA template and have less chance to frameshift or detach from 
the mRNA before they reach the stop codon. On the contrary, the plots in Figs. 5(e,f) indicate that, both /3g// 
and ratio r decrease with the rate ks of frameshifting. Considering large rate kg, ribosomes will have less chance to 
complete their translation processes correctly, but have more chances to frameshift. Figs. 6(b,c) show that, /3e// 
and ratio r also decrease with ribosome attachment rate uja- With the increase in attachment rate Wg, the mRNA 
template is bound by more ribosomes with incorrect peptide chains, and then the translocation speed of ribosome 
with correct peptide chain is reduced. Thus, ribosomes with correct peptide chains translocate along the mRNA with 
lower speed and then have more chances to detach from the mRNA or frameshift before they reach the stop codon. 
The plots in Fig. 6(e) show that /3g// may not change monotonically with the ribosome detachment rate ujd- When 
the detachment rate LOd is small, the increase in the rate uJd may help decrease the overall ribosome density along the 


mRNA, and consequently, increase the translocation speed of ribosomes with correct peptide chain. For large values 
of ^JJd^ the detachment rate of ribosome with correct peptide chain may increase rapidly, and then the synthesizing 
rate of correct peptide chain will decrease. This phenomenon implies that, one optimal ribosome detachment rate ujd 
may exist, with which the synthesizing rate of correct peptide chain reaches its maximum. 

The plots in Figs. 5(d), Fig. 6(a,d) show that the effective initiation rate cteff increases with the rate kg 
of frameshifting and the rate ujd of detachment, but decreases with the attachment rate Wa- With large detach¬ 
ment/attachment rate, the start mRNA codon has more chances to be unoccupied/occupied, therefore cteff is en¬ 
larged/reduced correspondingly. Fig. 5(d) also shows that with large values of frameshifting rate fcg, the effective 
initiation rate cteff tends to approach one limit constant. With large frameshifting rate feg, the entrance rate of 
ribosome to mRNA is only determined by the detachment rate uJd, the attachment rate uja, and the termination rate 
p. The overall density of the ribosome along the mRNA, with either correct peptide chain or incorrect peptide chain, 
is completely determined by uja,ujd,a, and /3, see [23-25]. The plots in Fig. 6 show that, although the translation 
effectiveness r always decreases with attachment rate Wa and detachment for some special cases the effective 
(correct) translation rate Peff may increase by increasing the ribosome detachment rate see Fig. 6(e). Therefore, 
detachment maybe one possible mechanism used by cells to increase the synthesizing rate of needed peptide chains. 


IV. CONCLUSIONS 

Recent experiments found that, during translation, ribosome may frameshift along the mRNA, and the reasons of 
frameshifting have been experimentally studied thoroughly [8-15, 17-20, 22, 40]. However, so far, no theoretical model 
has been designed to describe the gene translation process with ribosome frameshifting, especially the influence of 
frameshifting on translation, such as the effective translation rate and the efficiency (or effectiveness) of translation. 
In this study, a modified TASEP model is presented to describe the translation process with possible ribosome 
frameshifting. Similar to the usual translation process, i.e., translation without ribosome frameshifting, the probability 
density of ribosome along the mRNA template may have shock wave (or domain wall) and boundary layers. The 
effective translation rate, i.e., the rate of synthesizing correct peptide chain, increases with the translation termination 
rate /3, and decreases with the rate of frameshifting /cg, but may not change monotonically with the translation 
initiation rate a. Meanwhile, the translation effectiveness, which is defined as the ratio of effective synthesizing rate of 
correct peptide chain to the effective initiation rate of translation, increases with the termination rate, but decreases 
with initiation rate and the frameshifting rate of ribosome. At the same time, the influences of attachment/detachment 
of ribosomes to/from the main body of mRNA were also discussed. The results will be beneficial for the understandings 
of actual translation processes in cells, and the model presented may also be useful for the prediction of translation 
rate with reasonable chosen parameters. 
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TABLE I: Parameter values used in Fig. 3. In all calculations, N = 150, M = 50 and ks ~ 1 are used. 
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TABLE II: Parameter values used in Figs. 4, 5, and 6. Other parameter values used in the calculations are N = 150, M = 25 
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FIG. 1: (a) Normal translation process, in which the ribosome is always bound in the 0 frame of the mRNA. With one forward 
hopping of ribosome, a new amino acid carried by tRNA is added to the tail of peptide chain, (b) Incorrect translation, in 
which ribosome is bound in the —1 frame of mRNA. Compared with the normal ribosome in (a), the ribosome in (b) shifts 
one nucleic acid upstream along the mRNA (see the long vertical lines as reference position.) 
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FIG. 2: Modified TASEP to describe translation with ribosome frameshifting. The initiation of translation begins with the 
entrance of ribosome to the start codon, which corresponds to the left boundary of the lattice with site index 0. In any mRNA 
codon i (lattice site i), ribosomes may have four possible states, pi in (a) is the probability of the existence of a ribosome at 
codon i with a correctly synthesized peptide chain, while qf" are the probabilities of the presence of a ribosome at codon i of 
mRNA in ±1 frame. Given that ribosomes may bind to any mRNA codon, the peptide chain produced by a ribosome in 0 
frame may not be correct. For example, the ribosome at codon i may begin its translation from codon 1 < j < i- Thus, the 
peptide chain synthesized by it is also incorrect. The probability of the existence of such a ribosome at codon i is denoted 
by qi, see (b). Ribosomes in the 0 frame will leave mRNA from the stop codon [the codon N at the right boundary of the 
middle lattice line in (a)]. Generally, ribosomes in ±1 frame will not stop their translation processes until they reach the 3’ 
end of mRNA. M denotes the number of mRNA codons between the stop codon and 3’ end of mRNA. In the figures, a is 
the initiation rate of translation, /? is the termination rate of translation, and uJa and uJd are the attachment and detachment 
rates of a ribosome to and from mRNA, respectively. The elongation rate is denoted by fcs, which is assumed to be the same 
for ribosomes in any mRNA frame and at any codon. Finally, the rate of frameshifting, to either +1 frame or —1 frame, is 
assumed to be the same and denoted by kg. 
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FIG. 3: Typical examples of probabilities and their summation Pi — pi + Qi + + qf of ribosome along one mRNA 

with length N + M = 150 + 50 (see Fig. 2 for meanings of N, M and the corresponding probabilities.) The parameter values 
used in these plots are listed in Tab. I. Roughly speaking, ribosome translocation along the mRNA with frameshifting consists 
of two TASEPs, which are jointed at lattice site N. In (a,b,c), uja = u)d = 0, i.e., except at the start codon and stop codon, 
attachment and detachment of ribosome are not allowed. The legend given in (c) is right for all these nine figures. 
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FIG. 4: The effective translation initiation rate cteff := a(l — Po — go — — Qo) (a,d), effective (or correct) translation 

termination rate fieff '■= Ppn (b,e), and their ratio r := p^ff/cieff (c,f), as functions of the translation initiation rate a. In 
each figure, four typical examples are plotted, which are drawn by ‘x’, and ‘o’, respectively. I calculations, the initiation 
rate a changes from 0.01 to 0.5 with an increment of 0.01. The values of other parameters are listed in Table II. The main 
difference between (a,b,c) and (d,e,f) is that no frameshifting is allowed in (a,b,c), i.e. ks = 0, while in (d,e,f) ribosomes 
have nonzero frameshifting rate, i.e., ks > 0. 
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FIG. 5: The effective translation initiation rate Qe// (ajd), effective (or correct) translation termination rate /3e// (b,e), 
and their ratio r := (c,f), as functions of the translation termination rate (a,b,c), and the rate of ribosome 

frameshifting ks (d,e,f), where changes from 0.01 to 0.5 with an increment of 0.01, and ks changes from 0 to 0.01 with an 
increment of 0.002. See Table II for other parameter values. 
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FIG. 6: The effective translation initiation rate Oe// (ajd), the effective (or correct) translation termination rate /3e ff (b,e), 
and their ratio r := fi^ff/a^ff (c,f), as functions of the ribosome attachment rate LVa (a,b,c), and the ribosome detachment 
rate oJd (d,e,f), where cVa and ujd change from 0 to 0.01 with an increment of 0.002. See Table II for other parameter values. 


















